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We propose a simple and direct method for generating travelling wave solutions for nonlinear inte- 
grable equations. We illustrate how nontrivial solutions for the KdV, the mKdV and the Boussinesq 
equations can be obtained from simple solutions of linear equations. We describe how using this 
method, a soliton solution of the KdV equation can yield soliton solutions for the mKdV as well as 

' the Boussinesq equations. Similarly, starting with cnoidal solutions of the KdV equation, we can 

obtain the corresponding solutions for the mKdV as well as the Boussinesq equations. Simple solu- 
tions of linear equations can also lead to cnoidal solutions of nonlinear systems. Finally, we propose 
and solve some new families of KdV equations and show how soliton solutions are also obtained for 

| the higher order equations of the KdV hierarchy using this method. 

o" 

I. INTRODUCTION 

' Nonlinear integrable systems have been studied vigorously from a variety of points of view. Since the KdV equation 
was proposed the subject of integrable models has developed in many directions - the Lax representation, inverse 
scattering method, Backlund transformations, Zakharov-Shabat formulation, quantum integrability and quantum 
groups and so on 0, H, 0, H, d, @, 111 ■ Each of these developments is significant in its own right and has opened up 
new areas of research - see, e.g., [H, i, i, i, S S 0, 1, 0, El, El El, ES El, El, El, El El ll ■ The solution of a 
nonlinear integrable equation is normally obtained through the inverse scattering method jldl . 111 . E3] ■ I n some sense, 
this can be thought of as the "Fourier transform" method for nonlinear systems. Here one tries to associate with a 
given nonlinear integrable equation a pair of linear equations (with the dynamical variable of the nonlinear system 
£f} , as the potential) for which the "time" evolution of various quantities of interest is much simpler. By studying the 
"time" evolution of the scattering data in this linear system, one reconstructs the time evolution of the dynamical 
variable (the potential) through the method of inverse scattering. This is a powerful, indirect method for solving 

■ nonlinear integrable systems. There also exists a more direct way of solving a nonlinear integrable system that goes 
under the name of Hirota's bilinear method [f3 |. Here one introduces a non-conventional derivative (a derivative 

■ that acts forwards and backwards with a relative sign which is quite common in the study of quantum field theory) 
and rewrites the nonlinear equation as a linear equation in two variables which then leads to a direct solution of the 

i—i . original system. While both these methods are very powerful, they are quite technical. 

In this paper, we would like to propose an alternative direct method for solving nonlinear integrable systems that is 
, quite simple. In this method, one derives a nontrivial solution of a nonlinear integrable system from a known solution 
of a linear equation or another nonlinear equation. In the generic sense, therefore, this method can be thought of 

r% \ as a "Backlund" transformation method. However, unlike the conventional Backlund transformation method where 
one has to deal with zero curvature conditions, here the construction of the map is much more straight forward and, 
therefore, we believe that it is quite useful. This transformation method was already studied in connection with 
relativistic scalar field theories to generate nontrivial classical solutions to complicated theories from known ones in a 
simpler theory [2(|. In that context, it was also termed the "deformation" method, and in [2l| it was shown to work 
nicely for a variety of situations. However, the power of the method has led us to propose it as a useful direct method 
for obtaining solutions to nonlinear integrable systems. 

The paper is organized as follows. In section II, we illustrate the method within the context of a general Hamiltonian 
equation which is third order in space derivatives. In section III we show examples of how nontrivial travelling solutions 
for the KdV, the mKdV [1] as well as the Boussinesq [22| equations can be obtained from simple solutions of a third 
order linear equation. In section IV, soliton solutions for the mKdV as well as the Boussinesq equations are obtained 
from the soliton solution of the KdV equation. In section V, cnoidal solutions for the mKdV as well as the Boussinesq 
equations are derived from that of the KdV equation. We indicate briefly in this section how cnoidal solutions of 
nonlinear equations can be obtained from simple solutions of linear equations. The method proposed is then used to 
investigate new equations and the corresponding travelling wave solutions in section VI, and we conclude this work 
with some open questions in section VII. 
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II. ESSENTIALS OF THE METHOD 



Let us consider a general Hamiltonian equation in 1 + 1 dimensions of third order in the space derivatives of the 
form 

u t + (f(u)) x +au xxx = 0, (1) 

where the subscripts denote derivatives with respect to the particular variables and a is a constant. Here f(u) is a 
monomial in the dynamical variable u and the equation is linear/ nonlinear depending on the nature of this monomial 
function. Let us assume that this equation possesses a travelling solution of the form 

u T = u(kx — ujt), (2) 

where, in general, uj = uj{k) with the explicit form for the dispersion determined from the structure of the equation 
and in the case of a nonlinear equation, even the amplitude of the solution (height of the wave) is a function of k and, 
therefore, of u>. In the case of the KdV equation with /(it) = 3u 2 , a = 1 (the coefficients in the KdV equation can be 
rescaled to any value through appropriate transformations) 

u t + 6uu x + u xxx = 0, (3) 

for example, the one soliton travelling solution has the form 

u T = A sech 2 (kx-ut), lu = 4fc 3 , A = 2k 2 = ■ ( 4 ) 

This can also be rewritten in the more familiar form 

u T = 2k 2 sech 2 k(x - Ak 2 t) = - sech 2 k(x - ct), (5) 

where we have identified the speed of the travelling solution as c = 4/c 2 . We note here that if we transform u — > —u, 
the KdV equation takes the form 

u t - 6uu x + u xxx = 0, (6) 
for which the soliton solution (O has the inverted form 

u T = -2k 2 sech 2 k{x - Ak 2 t). (7) 
A travelling solution of (p} of the form (JSJ) then satisfies the equation 

- ^ + (f(u)) x + au xxx = 0, (8) 

which can be integrated to give 

u xx = ^- u-- f(u)+p = V'(u), (9) 
ak a 

where (3 is a constant of integration whose value can be determined from the asymptotic (boundary) conditions 
satisfied by the solution and we have introduced the notation 

VW - do) 

Here a constant of integration may be included in the definition of V{u), which in field theory is called a potential 

M- 

Let us next suppose that we have another dynamical system of third order with the dynamical variable v in 1 + 1 
dimensions and that we are interested in determining a travelling wave solution of the form 

v T = v(kx — ujt), (11) 
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where, in general, to = ui(k) will be a different function of k. As in the earlier case (see (O and (|10j0 . the travelling 
solution brings the dynamical equation for the system to the form 

v xx = V'(v) = ^1. (12) 
ov 

We now make the important observation that if there is an invertiblc map 

u = g(v), (13) 

such that HI 



V(9(v)) u v dg(v) 



n«) = fm *» = ^ (i4) 



then, Eq. ^ can be compared with Q12p to show that if a travelling solution for the first dynamical system is known, 
we can obtain a travelling wave solution for the second system simply as 

v T = v(kx — uit) = g~ 1 (u(kx — tot)). (15) 

It is important to note here that the relation (|13p is not required to map the first dynamical equation into the second 
(whose solution we are interested in constructing). Rather, it relates the derived potentials of the two systems. In 
this sense, it is different from a "Miura" map. 

We will show, through many examples, that this is a very powerful and useful method for constructing travelling 
wave solutions for nonlinear systems. As our discussion shows, the dynamical systems need not be integrable for the 
method to work. However, the utility of the method is clearly more significant for obtaining travelling solutions for 
integrable systems. In the next few sections, we will discuss in detail how such solutions are derived in a variety of 
cases. 

III. SOLUTIONS OF NONLINEAR EQUATIONS FROM LINEAR SYSTEMS 

In this section, we illustrate how this procedure works in generating nontrivial solutions for nonlinear equations 
starting with those for linear systems. Let us assume that f(u) = 0, a = 1 so that takes the form 

u t + u xxx = 0, (16) 

which is a third order linear equation and has a simple travelling solution of the form 

u T = A cos(kx + ut), lu = k 3 , (17) 

and the amplitude A is an arbitrary constant. In this case, ([9]) yields 

u xx = u + P = -k 2 u + (3 = V'(u). (18) 
k 

Comparing with the solution (|17|) . it follows that the constant of integration vanishes, (3 — 0, so that we have 

V(u) = ~ u 2 +1 , (19) 

where 7 is the constant of integration in the potential. 

If we are trying to obtain a nontrivial solution of say the KdV equation Q 

v t + 6vv x + v xxx = 0, (20) 

of the form 

u T = v{kx — u>t), (21) 
we can follow the discussion in (j9j)- (fT0|) to construct 

V(v) = ^v 2 -v 3 +/3v + j, (22) 
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where 7 is the integration constant in the potential. 

Next, let us note that if we choose the constant of integration 7 = A 2 k 2 /2 for simplicity, then under the map 



u(kx — ujt) = g(v(kx — ujt)) = A cos ^sech 1 f^J ^g^j^j > (^3) 



we obtain 



^1 = 2fcV - , 3 . (24) 

{g'{v)f 

This can be compared with (f2"2"|) and we note that we can identify V(v) as in (|14[) provided 

/3 = 7 = 0, £ = 4fc 3 . (25) 

In this case, the solution for the KdV equation can be obtained as (see (jT5j) ) 

v T = g^ 1 {u{kx - uit)) = 2k 2 sech 2 k(x - Ak 2 t), (26) 

which we recognize as the one soliton solution of the KdV equation ([5]) and we note here that the vanishing of the 
constants of integration in (|25p is a reflection of the fact that the soliton solution and its derivative vanish asymp- 
totically. In this case, we have constructed this nontrivial solution of the KdV equation from a simple trigonometric 
solution of the linear equation p6[) . 

On the other hand, if we want to construct a solution of the mKdV equation with the negative sign 

v t - 6v 2 v x + v xxx = 0, (27) 

of the form 

v T = v(kx — uit), (28) 
then, we can proceed as in ©-(US]) an d determine 

V{v) = ^v 2 + l -v A + (lv + ^. (29) 
Here /?, 7 are the two constants of integration. It is clear that if we choose as before 7 = A 2 k 2 /2, then under the map 

u{kx — cut) = g(v(kx — uit)) = ±Acos ^tanh -1 yrjj 5 (30) 

we obtain 

(g'(v)f 2 2' 
which can be compared with (|29|) . We note that we can identify V(v) as in (fM]) provided 

/3 = 0, 7 = y, u = -2k\ (32) 

In this case, the solution of the negative mKdV equation (|27p is obtained to be 

v T = g~ x (u(kx - uit)) = ±fctanhfc(a; + 2k 2 t). (33) 

This is the well known soliton solution of the negative mKdV that is not localized, but here we have derived it from a 
simple trigonometric solution of a linear equation. (Parenthetically we note that a nontrivial constant of integration 
7 in this case signals that the solution does not vanish asymptotically (since it is not localized), but the derivative 
does.) 

The travelling solution for the (positive) mKdV equation 

v t + 6v 2 v x + v xxx = 0, (34) 
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can also be constructed in a similar manner as follows. For a travelling solution of this equation of the form 

v T = v(kx — u)t), (35) 

following the steps ([9|)- (fT0|) . we obtain from |34|) 

V{v) = ^v 2 ~ l -v A +{3v + l- (36) 

Here (3 and 7 are the constants of integration and this potential can be compared with that in (f!?9"]) . We note that, as 
in the earlier cases if we choose 7 = A 2 k 2 /2, then under the map 

u(kx — tot) = g(v(kx — tot)) = ±Acos ^sech -1 (rJJ 1 (37) 



we obtain 



V(g(v)) k 2 v 2 1 



{g'{v)f 2 2 

which can be compared with (f3"6"|) . It is clear that we can identify V(v) as in (|14[) provided 

/? = 7 = 0, w = /c 3 . (39) 
In this case, we can determine the travelling solution for the mKdV equation to be 

v T = g^ 1 (u(kx - cot)) = ±k sech k(x - k 2 t), (40) 



7T v\ (38) 



which represents the one soliton solution of the mKdV equation. We note here that the two mKdV equations, (|27|) 
and (|34| . arc invariant under v — > — w, which is the reason for the solutions with both signs in these two cases. 

Finally, let us note that we can also construct a travelling solution for the Boussinesq equation in a very simple 
manner from the solution (|17[) as follows. First of all, we note that the Boussinesq equation is a higher order equation 
of the form 

v tt -v xx ~3 (v 2 ) xx - v xxxx = 0. (41) 
Nonetheless, if we assume that it has a travelling solution of the form 

v T = v(kx — cot), (42) 
then, following the steps in (|9j)- (fT0|) . it can be brought to the form (fT2|) with 

^) = ^((l) -ljv 2 -v 3 + /3v + j. (43) 



Furthermore, as we have already seen, under the map (|23D . 

V(g(v)) 



(g'(v)f 



= 2k 2 v 2 - v 3 , (44) 



which can be compared with (|4"3")) and it is clear that in this case we can identify V(v) as in (TH| provided 

/3 = 7 = 0, lo = ±fei/l + 4fc 2 . (45) 
The nontrivial travelling solution for the Boussinesq equation now follows from (|23|) to be 



v T = g~ l {u{kx - Lot)) = 2k 2 sech 2 k(xT Vl + 4fc 2 t). (46) 

This is indeed the single soliton solution of the Boussinesq equation and we note that unlike the soliton solution for 
the KdV (or the mKdV) equation, that of the Boussinesq equation is bi-directional. 
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IV. SOLITON SOLUTIONS FROM SOLITONS 



In the last section, we tried to illustrate the details of the method of our proposal by constructing nontrivial soliton 
solutions of nonlinear integrable systems from a simple trigonometric solution of a linear equation. In this section, 
we continue the exposition of the method by showing how one can construct soliton solutions for other integrable 
systems starting with a soliton solution of a given nonlinear integrable system. For simplicity, let us assume that we 
start with the one soliton solution (|26[) (or (J3J)) 

u T = 2k 2 sech 2 k(x - 4k 2 t), (47) 

of the KdV equation j3]) 

Ut + 6uu x + u xxx = 0. (48) 

The potential for the KdV equation following the procedure (|9j)-(fT0|) has already been constructed in (|22|) . However, 
for the soliton solution in (|4"7| , we note that both the solution as well as its derivatives vanish asymptotically which 
determines the constants of integration to be trivial (which we have already noted in the last section). As a result, 
for this case, we have 

V(u) = 2k 2 u 2 - u 3 . (49) 



Let us construct a travelling solution of the negative mKdV equation (|27j) starting from the solution (|47|) . We have 
already seen in (f2"5)) that for a travelling solution of (J27J) 

V(v) = v 2 + - v A + (3v + 7. (50) 

It is clear now that under the map 

u(kx - ut) = g(v(kx - ut)) = 2(k 2 - v 2 ), (51) 

we can identify 

= fcV + i , 4 = V(v), (52) 

(g'(v)) 2 2 2 

provided we have 

/3 = 0, l=\k\ u = -2k\ (53) 

This, in turn, leads to the soliton solution 

v T = g -1 (u(/ca; - ujt)) = ±fctanhfc(a; + 2k 2 t), (54) 



which we have derived earlier in 

For the (positive) mKdV equation ([34]) . on the other hand, we have already determined (see (|36|) ) 



We note that under the map 



we can identify 



provided we have 



V(v) = ^-v 2 ~^v 4 + pv + j. (55) 



u(kx — ujt) = g(y{kx — ut)) = 2v 2 , (56) 
!^ = i JbV - i «* = V(f,), (57) 
/3 = 7 = 0, u = k 3 . (58) 
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The inverse of the map (|56j) now leads to the solution 

v T = g~ 1 (u(kx — u)t)) = ±k sech k(x — k 2 t), (59) 



which we have also seen in (|40j). 

At this point, it is worth comparing our construction with the conventional Backlund transformation for KdV. We 
note that the Miura transformation 

u = (3v 2 ± y/^Pv x , (60) 

maps the KdV equation ([3]) to the positive as well as the negative mKdV equations for /3 = 1 and j3 = — 1 respectively. 
For (3 = 1, the transformation is clearly complex and it is worth noting that the Miura transformation alone cannot 
determine the time dependence of the solution (dispersion). For that one has to combine the Miura transformation 
(|60|) with the evolution equation for mKdV to obtain the Backlund transformation for KdV, which has the form 

v x = ±_ ^=^ {u- (3v 2 ) , 

v t = T-^{u xx ±2y^(3u x v + 2u 2 -2(3uv 2 y = ±1. (61) 

The compatibility of the two equations leads to the KdV equation and mapping solutions between the two systems 
corresponds to working with the above set of equations which is rather nontrivial compared to our proposal. Fur- 
thermore, it seems reasonable that Backlund transformations will map corresponding classes of solutions of the two 
systems into one another (namely, solitons to solitons), but it is not a priori obvious whether they can map between 
classes of solutions (and even if it can be, it seems rather nontrivial). Here our proposal seems to be more versatile 
in that it can map solutions between different classes. In this section, we have shown how soliton solutions of KdV 
can generate soliton solutions of mKdV. In the next section, we will indicate how soliton solutions of KdV can also 
generate cnoidal solutions of mKdV. 

Let us conclude this section by indicating how the soliton solution for the Boussinesq equation (|41[) is obtained 
from (|4"T)) . We have already seen in that a travelling solution of the Boussinesq equation leads to 

v ( v ) = tH ( T ) - 1 I v 2 - y3 + P V + 1- (62) 

It now follows that under a map 




we can identify 



provided we have 



u(kx — tot) = g(v(kx — tot)) = v, (63) 
V(g(v)) 



(g'(v)f 



2/sV - v d = V(v), (64) 



/? = 7 = 0, w = ±k^l + Ak 2 . (65) 
In this case, the inverse map leads to the one soliton solution of the Boussinesq equation, namely, 



v T = g- 1 (u{kx - Zjt)) = 2k 2 sech 2 k(x =p Vl + 4fc 2 t), (66) 
which is what we have derived earlier also. 



V. CNOIDAL SOLUTIONS 



Nonlinear integrable equations have a rich variety of nontrivial travelling wave solutions. In the last section we 
discussed one of these types of solutions, namely, the soliton solutions. There are also the cnoidal solutions which are 
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more general than the soliton solutions and reduce to the soliton solutions in particular limits. Thus, for example, for 
the KdV equation ^ 

u t + 6uu x + u xxx = 0, (67) 

it is known that a cnoidal solution has the form 

u T = 2mk 2 cn 2 (kx-Lot,m)+B, lu = k(6B + 4(2m - l)fc 2 ), (68) 

where B is an arbitrary constant and to is a constant lying between < m < 1. The cn function is an elliptic 
function. There are other such elliptic functions which are also solutions of the KdV equation and they satisfy 
interesting properties. In particular, when m = 1, the cn function reduces to the "sech" function and if we choose 
m = 1, B = 0, then the cnoidal solution of KdV in (f6"5|) reduces to the one soliton solution ([5]). 

Unlike the soliton solutions of the KdV equation, neither the cnoidal solutions nor their derivatives vanish asymp- 
totically. Therefore, for such a travelling solution, following (j9|)- (TlO|) . we can determine (see also (|22jl ) 

V(u) = ^- u 2 -u 3 + (iu + ll (69) 

where lo is defined in (|68p and /3, 7 are constants of integration whose values can be determined from the properties 
of the cnoidal solutions. In particular, it is known that the "potential" in (|69p can be factored in terms of the three 
roots of the equation as 

V(u) = (B + 2(m - l)fc 2 — u)(B — u) (B + 2mk 2 - u) . (70) 



Given the cnoidal solution of the KdV equation in (|68j) . let us now construct the cnoidal solution for the mKdV 
equation (|34|) following our procedure. First, we recall from (|36j) that for a travelling solution in this case, we can 
write 

V{v) = ^v 2 - l -v A + Pv + l- (71) 

Furthermore, under a map of the form 

u(kx - ut) = g(v(kx - ut)) = 2v 2 + B, (72) 



YMO^ = _i m r m _ i) k i + i ( 2m - l)k 2 v 2 --v 4 = V(v), (73) 



(3 = 0, 7 = ~ m(m- l)fc 4 , i = (2m-l))t 3 . (74) 



we can identify 



provided we have 



The cnoidal solution for the mKdV equation can now be obtained from the inverse of the map (|72|) 

v T = g^ 1 (u(kx — ut)) = ±^/m k cn(kx — uit, m). (75) 

where lu is given by (|74p . 

For m = 1, these reduce to the soliton solutions of the (positive) mKdV equation in (|4"D")) . Note also that the 
constant introduced in the map (f72|) corresponds to one of the roots of the "potential" in (|70|) . We could instead have 
defined the map with any other root of the potential to obtain a different cnoidal solution of the mKdV equation. 
In particular, the other roots in ([70)) will lead to solutions which will reduce to the soliton solution of the (negative) 
mKdV in (|33|) for m = 1. As we have noted in the previous section, our method can lead easily from one class of 
solutions to another. So, for example, starting from the one soliton solution of KdV (see (|47|) and (1491 ). we note that 
the map 



uikx — tot) = g(v(kx — Lot)) = 2fc 2 sech 2 



cn 1 ( ± , to 

mk 



(76) 



leads to (fTB")) . The inverse map then generates the cnoidal solution ([75")) starting from a soliton solution of KdV. This 
demonstrates the versatility of our method. 
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The cnoidal solution of the Boussinesq equation (|4Tj) can also be constructed following our procedure. For a 
travelling solution of the Boussinesq equation, we have seen in (|43[) that we can write 



V(v) 



- 1 }v 2 - v 3 + f3v + 7. 



In this case, we see that if we define the map 

u(kx — lot) = g(v(kx — lot)) = v, 

then we can identify 

V(g(v)) 



(g'(v)Y 



(B + 2(m - l)k z —v)(B — v) (B + 2mk 2 



V(v), 



(77) 



(78) 



(79) 



provided we have 



4m(l - m)k 4 + 4B(1 - 2m)fc 2 - 3B 2 

4Bm(m - l)fc 4 + 2B 2 (2m - l)fc 2 + B 3 , w = ±fci/l + 6-B + 4(2m - l)fc 2 . 



/3 
7 

The inverse of the map (|78[) now yields the cnoidal solution of the Boussinesq equation 

v T = g~ 1 (u(kx — tot)) = 2mk 2 an 2 (kx — u)t 1 m) + B, 



(80) 



(81) 

where u> is given by (|80p . 

Finally, let us close this section by noting that just as the soliton solutions of nonlinear systems can be obtained 
from travelling solutions of a linear system such as (|16p . so can the cnoidal solutions be obtained from travelling 
solutions of a linear system. Without going into details, let us indicate here how this can be carried out in obtaining 
the cnoidal solution of the KdV equation ((3J starting with the simple trigonometric solution, (fT7|) . of the third order 
linear equation P^)l . Let us choose 7 = A 2 k 2 /2 in p^|) (as before) and note that under the map 



i(kx — uot) = g(v(kx — Lot)) = ^4cos 



v — B 
2mk 2 ' 



(82) 



we can make the identification (|T9"|) with V(v) given in ([^)) . provided we choose the constants (3 and 7 as in (fBT))) . and 

lo = fc(6B + 4(2m- l)fc 2 ). (83) 



In this case, the inverse of the map (f82|) leads to the cnoidal solution (|68|) of the KdV equation. Similarly, the cnoidal 
solution (|75|) for the mKdV equation (|34p can be generated from this solution of the linear equation under the map 



u(kx — Lot) = g(v(kx — Lot)) = A cos 



m k 



, m 



(84) 



This demonstrates that cnoidal solutions, like the soliton solutions, can also be obtained from simple solutions of 
linear systems, using our method. 



VI. SOME NEW NONLINEAR EQUATIONS AND THEIR TRAVELLING SOLUTIONS 

So far we have used our method to construct (known) solutions of known equations. We now turn our attention to 
applying this method to solve new nonlinear equations. We first consider a dynamical system described by 

v t + o(l +p)(l + 2p)v 1/p v x + v xxx = 0, (85) 

where a and p are arbitrary, real constants. We name this the pKdV equation, and note that for a = p = 1 it reduces 
to the KdV equation while for a = ±2,p — 1/2, it leads to the two mKdV equations. We analyze this new equation 
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with the method discussed in section II and construct its travelling wave solutions starting from that of the KdV 
equation. For a travelling solution of (|85[) . we can determine (following (|9"))- (flT))) ) 

V(v) = ^-v 2 ~ ap 2 v (1+2p}/p + 0v + 7. (86) 
2k 

We note that under the map 

u(kx — Lot) = g(v(kx — uit)) = av 1 ^, (87) 
with the form of u given in ((4]). we can identify 

YM^l = 2k 2 p 2 v 2 - ap 2 v^ 1+2p ^ p = V(v), (88) 
O'O)) 

provided we have 

= 7 = 0, Z) = 4&V- (89) 
In this case, the inverse of the map (|87|) leads to the solution 

w T = 5 " 1 (u(fca; - St)) = I — J sech 2p k(x - 4k 2 p 2 t). (90) 

We note that a real solution exists for any value of p when a > 0. However, if a < 0, the values of p for which real 
solutions exist are restricted. For the particular value a = l,p = 1 as well as for a = 2,p = 1/2, this solution nicely 
reduces to (|26p and (jlO")) respectively. 

Let us next construct the solutions of a different family of nonlinear equations as another illustration. We start 
with the negative mKdV equation (|27p and we consider the kink-like soliton solution u T = ±fctanhfc(x + 2k 2 t). We 
are then led to 

V{u) = \{u 2 -k 2 ) 2 . (91) 

We introduce the map 

u(kx — w m ,a t) = g(v(kx — w m ,a t)) = iik cos (a arccos(w) — mir) , (92) 
where a is a real parameter, and m takes integer or semi-integer values. For m semi-integer we can identify 

Vi9{V)) - lk2 ;(l-v*)T 2 M = V 1>a (v), (93) 



and for m integer we have 



(g'(v)) 2 2 a 2 
V{9(v)) Ik 2 2 2 ~ 

(94) 

where we have denoted T a (v) = cos(a arccos(u)) and U a (v) = sin(a arccos(?;)), for simplicity. For a integer, T a (v) and 
U\+ a {v) are the Chebyshev polynomials of the first and the second kinds, respectively. We apply our procedure to 
the two families of equations 

ut + (/<,«.(«))' v x +v xxx =0, (95) 

where fi. a (v) is a function of v, and i — 1 or 2 correspond to the two possibilities in (|93|) and (|94| respectively. We 
name this the KdV(i,a) family, and note that for travelling waves we have 

y^/,>)+i = (iW)', (96) 



with solutions 

v T = g~ 1 (u(kx — Wi a t)) = cos [ — (6(kx — Wi a t) + mir) I , (97) 

\a J 

where the velocitiy uji_ a depends on m being semi- integer (i = 1) or integer (i = 2). Here 9{y) is the principal 
determination of arccos (tanh(y)) , and m now identifies one among the several solutions the equation may comprise: 
for m semi-integer we have m = 1/2, 2a — 1/2, and for m integer, m = 0, 1, 2a — 1, for a being integer or semi- 
integer. In the figures below one identifies these possibilities very cleary. To make these illustration more explicit, let 
us now examine in detail the different families that we have for a integer or semi-integer. 
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A. KdV family for a semi-integer 

In order to obtain the explicit form of (/j, a (w))' and the related Wi <a , let us write Vi, a (v) as a function of its zeroes. 
For a semi-integer we have 

V ha (v) = «)(1 - v 2 ) J] (« + 4 .a) 2 , (98) 

and V2, a {v) = Vi,a( — w). In the above expression we have set 

Z la = COS ( ^-=^7T ) . (99) 



Using J96J and JSSJ) we obtain 



2a 



(/!»)' = (100) 



where we have used «?! „ = —k 3 /2a 2 . The polynomial functions (f9"8")l provide a family of KdV equations which 
presents one sech-likc soliton solution and N = a — 1/2 different pairs of tanh-like soliton solutions, for each member 
of the family of equations. We illustrate some of the many possibilities with the cases a = 1/2, 3/2, and 5/2 : 

(/i,i/a(«))' = 6 k 2 v, w hl/2 = -2k 3 , (101a) 




FIG. 1: Plots of Vi, (t>) (top panel) for k = 1, and a = 1/2 (left panel), a = 3/2 (middle panel), and a — 5/2 (right panel), and 
the corresponding sech-like and tanh-like solutions (botton panel) for t = 0. 
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B. KdV family for a integer 

In the case of a integer, we have distinct families of KdV equations, one for i = 1 and the other for i = 2. For this 
reason, below we deal with the two possibilities separately. 

1. The case i = 1 

In order to obtain the explicit form of (fi : a(v))' and the related w\^ a for a integer, let us write V\ ta {v) as a function 
of its zeroes. We have two possibilities: 

(i) a even 

fc» = _ ^ {z{ ff. (102) 

3 = 1 

(ii) a odd 

^ Q («) = ^^ V 2 (1-, 2 ) J] (^-(4,J 2 ) 2 , (103) 



where 



Using (dU), (TUg) and CnnD we obtain 



(/i,W)' = y-(^W)", (105) 

with wi i£1 = A; 3 for a odd, and Wi <a — ~k 3 for a > 2 even, with |tux ja +2| < |wx, a |. The polynomial functions (|102|) 
and (|103p provide families of KdV equations which present two sech-like solutions and N = a — 1 pairs of tanh-likc 
solutions, for each member of the family. We illustrate some of the many possibilities for the cases a = 1, 2, 3 : 

(fi,x(v))' = 6k 2 v 2 , u5 1(1 = fc 3 , (106a) 

(f l 2 ( V )Y = -12 k 2 v 2 + 15 k 2 v 4 , wi a = -1-25 fc 3 , (106b) 

(A»)' - 22 fc 2 z, 2 - k * v * + 115 fc 2 ^ ~ i 3 = fc3 _ (1Q6c) 
In Fig. 2 wc plot Vi, a {v) and the related solutions, for some integer values of a. 

2. The case i = 2 

In order to obtain the explicit form of (f2,a(v))' and the related u>2,a, let us write V2, a (^) &s a function of its zeroes. 
Here we have: 

(i) a even 

= ^^, 2 JJ^ " (^,a) 2 ) 2 - (107) 

i=l 

(ii) a odd 

9 2a-3 J.2 f^ 1 )/ 2 

^2» = ^— ^ I] (- 2 -(4a) 2 ) 2 , (108) 
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FIG. 2: Plots of Vi, a {v) (top panel) for k = 1, and a — 1 (left), a = 2 (middle), and a — 3 (right), and the corresponding 
solutions (botton panel) for t = 0. 



where 




(109) 



Using (|96l), (fT07j) and (fT08|) we obtain 

(/2,.(«))' - (HO) 

and u)2,o = k 3 for a even, u5 2 ,i = —2 fc 3 , and u>2, a — — k 3 for a > 1, odd, with |w2,a+2| < 1^2, o|- The polynomial 
functions (| 1 07[) and (|108|) provide families of KdV equations which present N = a pair of tanh-like solutions, for 
each member of the family. In Fig. 3, we illustrate some of the many possibilities with the cases a = 1,2,3, for 
w 2m = ^2 fc 3 ,/c 3 ,— 1.11 fc 3 , respectively. 

VII. CONCLUDING REMARKS 

In this paper, we have proposed a simple and direct method for finding travelling solutions of nonlinear integrable 
systems. We have illustrated the method by explicitly constructing soliton solutions of the KdV, the mKdV and the 
Boussinesq equations starting with a simple trigonometric solution of a linear equation. We have also shown how 
the method can be used profitably to construct soliton solutions from a known soliton solution of a given integrable 
equation. In this connection, we have constructed the soliton solutions of the mKdV and the Boussinesq equations 
starting from the known soliton solution of the KdV equation and have compared our method to the conventional 
Miura/Backlund transformation methods. The method applies to other types of travelling solutions as well and we 
have emphasized this by constructing the cnoidal solutions of the mKdV and the Boussinesq equations starting with 
that of the KdV equation. We have also indicated how the cnoidal solutions can be obtained from simple solutions 
of linear systems. We have chosen these examples (of mostly known solutions) primarily to emphasize clearly how 
the method works in practice and how simple and useful it is. (One can contrast this with the technicalities that go 
into the construction of each of the known solutions.) There are many other examples that we have not been able to 
describe because of space limitation. 

The procedure is cleary very powerful and useful. As we have explicly shown, it works very nicely to contruct and 
solve distinct families of given equations. In particular, we have constructed soliton solutions for pKdV equations. 
We have also shown how to construct equations which support sech- and tanh-like solutions. As we know, this is in 
contrast with the standard situation, where the KdV equation supports scch-like soltions, and the mKdV, tanh-like 
solutions. The results of Sec. VI show that one can find sech- or tanh-like solutions for the same equation, depending 
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on the boundary condition one uses to get the solutions. We can also check that, for the KdV hierarchy (we simply 
quote the result here without going into the details), the one soliton solution for the n-th equation in the hierarchy 
can be constructed using our method and has the form 

4 2 " +1) = 2fc 2 sech 2 k(x - {2k) 2n t), (111) 

where for n — 1 (the KdV equation ([3])) the solution corresponds to the soliton solution 

The results of the present work compel us to ask whether we can study multi-soliton solutions of integrabe systems, 
and to solve higher order equations in a systematic manner. It would be interesting, for example, if our method can 
shed light on whether a soliton solution exists for an equation of the form 

ut + auu x + Suxxxxx = 0. (112) 
These are some of the questions that arc presently under study and we hope to return to them in a future publication. 
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